Double Beta Decay: Theory, Experiment, and 
Implications 
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1 Introduction 



Double beta decay is a rare spontaneous nuclear transition in which the nu- 
clear charge changes by two units while the mass number remains the same. 
It has been long recognized as a powerful tool to study lepton number con- 
^\ [ servation in general, and neutrino properties in particular. Since the lifetimes 

of double beta decays are so long, the experiments on double beta decay are 
very challenging and have led to the development of many generally valuable 
techniques to achieve extremely low backgrounds. 
£vq \ For the (3f3 decay to proceed, the initial nucleus must be less bound than 

f^) ■ the final one, but more bound than the intermediate nucleus. These conditions 

l/"") | are realized in nature for a number of even-even nuclei (and never for nuclei 

■ with an odd number of protons or neutrons). Since the lifetime of the (3(3 decay 

is always much longer that the age of the Universe, both the initial and final 
nuclei exist in nature (some of the actinides being the only exceptions). In 
many of the "candidates" the transition of two neutrons into two protons 
is energetically possible, with the largest Q value just above 4 MeV. In a 
O \ few cases the opposite transition, which decreases the nuclear charge, is also 

possible, but the Q values are typically smaller. 

The nuclear (3(3 transition can proceed in several ways. One of them, the 
2v decay 



(Z,A) (Z + 2,A) +e l +e 2 + D ei + D e2 (1) 

conserves the lepton number, while the other one, the Ov decay 

(Z,A)^(Z + 2,A) + e^ +e 2 - (2) 

violates lepton number conservation and is therefore forbidden in the stan- 
dard electroweak theory. The prospect of discovering this neutrinoless double 
beta decay mode is the driving force of most of the interest in this field. It is 
the possible window into physics "beyond the Standard Model" . 

Double beta decay has been and continues to be a popular topic since 
first discussed by Maria Goeppert-Meyer in nineteen thirties. There have 
been numerous earlier reviews, beginning with the "classics" by Primakoff 
and Rosen |j} , Haxton and Stephenson || and Doi, Kotani and Takasugi || , 
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to the more recent ones, often devoted to particular aspects of the (3(3 decay 
[|J|^,0|8[|Q. Many details are also described in the monograph |nj|. The 
Review of Particle Physics regularly summarizes the most recent exper- 
imental data. 

Double beta decay in all its modes is a second order weak semileptonic 
process, hence its lifetime, proportional to (Gf cos 0c)~ 4 , is so very long. 
(Here Gf — 1.166 x 10~ 5 GeV -2 is the Fermi coupling constant, and 9c is 
the Cabbibo angle.) The neutrinoless decay can be mediated by a variety of 
virtual particles, in particular by the exchange of light or heavy Majorana 
neutrinos. The decay amplitude then depends on the masses and coupling 
constants of these virtual particles. Independent of the actual mechanism 
of the Ov (3(3 decay, its observation would imply that neutrinos necessarily 
have a nonvanishing Majorana mass fl3|| . In fact, if the Oi/ decay is actually 
observed, and its rate measured, one can obtain, at least in principle, a lower 
limit on that mass . 

However, so far no Qv decay has been observed. This means (barring arti- 
ficial complete cancellation of the amplitudes which we dismiss as unnatural) 
that the upper limit of the decay rate can be interpreted as an indepen- 
dent limit for each of the possible amplitudes of the decay. In particular, we 
can obtain the limit on the properties of light and heavy virtual Majorana 
neutrinos. Below we concentrate on the decays mediated by these particles. 
(Other possibilities, e.g., the decays mediated by the new particles predicted 
by supersymmetry, are discussed in p|^0[.) 

The [3(3 decay with the Majoron emission, Ovx mode, 

(Z,A)^(Z + 2,A) + e^ + X . (3) 

belongs to the category of the lepton number violating decays, even though 
the lepton number is formally conserved when \ is assigned the lepton number 
—2. The hypothetical scalar particle x, which must be in this case light 
enough to be emitted in the (3[3 decay, is usually associated with spontaneous 
breaking of the B — L symmetry [|l5|Jl6[] . 

Empirically, it is easy to distinguish between the three decay modes listed 
above, provided the electron energies are measured. The electron sum energy 
spectra are determined by the phase space of the outgoing leptons and clearly 
characterize the decay mode, as schematically illustrated in Figure 0. (Geo- 
chemical or milking experiments, however, cannot distinguish between the 
different [3(3 modes as they determine only the total decay rate.) 

There are two distinct groups of theoretical issues associated with the 
interpretation of the (3(3 decay experiments. The particle physics issues deal 
with the expression of the decay rate in terms of the fundamental parameters, 
such as the neutrino masses and mixing angles, coupling constants in the weak 
interaction Hamiltonian, etc. This group of problems involves also the relation 
of the (3(3 decay to other processes, such as neutrino oscillations, direct mass 
measurements, and searches for other lepton number violating processes. 
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Fig. 1. Schematic sum electron spectra of the three (3(3 decay modes. Each is nor- 
malized arbitrarily and independently of the others. The abscissa is the ratio K/To 
of the sum-electron kinetic energy divided by its maximum value. 



The other, essentially decoupled, set of problems involves the nuclear 
structure issues associated with the (3(3 decay. The decay rate is expressed in 
terms of nuclear matrix elements (NME) which have to be evaluated. One 
would like to know, first of all, their value and its uncertainty. This area of 
research has attracted lots of attention, and there are many, often conflicting, 
evaluations available in the literature. Unfortunately, there is no simple way 
of judging the correctness and accuracy of the evaluations of the nuclear ma- 
trix elements for the neutrinoless decay. Comparison to the experimentally 
known rate of the 2v (3(3 decay rate is often invoked in that context as a test 
of the ability of the nuclear model to describe the related phenomena. It is 
not clear, however, if this is indeed a valid test. For example, if one assumes 
that the Qv decay is mediated by the exchange of a heavy particle (whether 
this exchanged particle is a heavy neutrino or not), the corresponding inter- 
nucleon potential is of short range, and additional issues involving nucleon 
structure, irrelevant for the 1v decay, play an important role. Another, less 
fundamental but in praxis perhaps more important example deals with the 
dependence of the NME on the number of single nucleon subshells included 
in the calculation. For the 1v decay, where only the Gamow- Teller operator 
err plays a role, it is clearly sufficient to include just the states within the 
valence oscillator shell. It is less clear that the same truncation is sufficient 
for the correct description of the Ov decay. 

The experimental study of the (3(3 decay presents a formidable challenge 
since the goal is to detect a process with a half-life in excess of 10 25 years 
(the present best limit for the Qv decay). The (3(3 decay must be detected in 
presence of an inevitable background of a similar energy caused by trace ra- 
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dioisotopes with half-lives 15 or more orders of magnitude shorter. Thus, the 
optimum separation of the signal from background, combined with the re- 
quirement of having kilogram quantities of the source isotopes, characterizes 
the present day experiments. 

The past and current experiments are still relatively modest in size, and 
therefore also in complexity and cost. Given the importance of the search 
for the neutrinoless decay, ambitious plans, involving much larger amounts 
of the source nuclei, are considered. Naturally, the larger source mass will 
be beneficial only if it is accompanied by the corresponding reduction of 
the background. The future projects will be therefore inevitably much more 
complex, and will involve larger groups of researchers. With them, the field 
of (3(3 decay, which competes already with the other experiments described 
in this book in importance, will also compete in size and cost. 



2 Lepton number violation 

With the usual assignment of the lepton number, L(l~) = L(v) = —L(l + ) = 
—L(p) = +1, Ov (3(3 decay represents a change in the global lepton number 
by two units, AL = 2. In that respect its observation would be related to the 
attempts to detect 9 e from the sun, or of v e from nuclear reactors. Both of 
these latter processes represent a kind of u v <-> v oscillations" , and also are 
possible only for massive Majorana neutrinos. 

For light Majorana neutrinos the lepton number conservation is irrelevant 
and the Ov decay is hindered only by the helicity mismatch. However, the 
"antineutrino" born in association with one of the e~ in the Qv decay is not 
fully righthanded, but has a lefthanded component of amplitude ~ m v /E v . 
This lefthanded piece can be absorbed by another neutron which is converted 
into a proton and the second e~ is emitted. Similar consideration would 
govern the above mentioned v <-» v oscillations. The word "oscillations" in 
this context is a misnomer, however, since the process (if it exists) would 
proceed without an oscillatory behavior JItJ . 

The expected branching ratio for the "wrong" neutrinos at low energies, 
relevant for the sun or nuclear reactors is [fl8| 

m 2 rr 9N 

(«) 

where the numerical factor was derived for m v ~ 1 eV, E v ~ 5 MeV, and 
the ratio of cross sections put to unity. Since the Qv (3(3 decay is presently 
sensitive to such neutrino masses, one cannot expect a signal for this kind of 
v «-> v oscillations until similar sensitivity is achieved, i.e., not anytime soon, 
if ever. 

However, it is also possible that v e from the sun are produced in a more 
complicated, but possibly more efficient way. Let us assume that a transition 
magnetic moment /i ej ; connects the lefthanded v e> L with a righthanded v^r of 
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a different flavor, which can subsequently oscillate (by the vacuum or matter 
enhanced oscillations) into the righthanded and thus observable v e ^ 1 i.e., 
when neutrinos propagate in a transverse solar magnetic field B± one or both 
of the sequences v e ^ —* H,r — > v & ,r or v e L — ► — * ^ei?, occurs. Such process 
requires that the magnetic conversion, which is possible only for the massive 
Majorana neutrinos and which depends on the product ^ e ,iB±, and the flavor 
oscillation, which depends on Am 2 and sin 2 2(2, are both present. There is no 
obvious relation between this process and the neutrinoless (3(3 decay, except 
that both require the existence of the neutrino Majorana mass term. (This 
brief discussion of the magnetic conversion is highly simplified. In reality, the 
transition magnetic moments ought to be written in terms of mass eigenstates 
&) 

Finally, tight experimental limits exist on the total lepton number violat- 
ing processes which involve both electrons and muons (see fl^| ), such as the 
muon conversion 

ti- + (Z,A)^(Z-2,A) + e + , (5) 

and the muonium-antimuonium conversion 

/i + e~ — > [iT e + . (6) 

The relation of these processes to the (3(3 decay is, however, not well estab- 
lished. 



3 Particle physics aspects 

In this section we shall consider how the rate of the neutrinoless (3(3 decay is 
related to the unknown parameters of the neutrino mass matrix and to the 
phcnomcnological parameters describing a generalized semilcptonic charged 
current weak interactions H\y'- 

Hw = ^[J£(M+ Q + ^4J + J%(vM+ a + \M+ a )] + H.c. , (7) 

where Jl(r) and Mum are the lepton and quark left(right)-handed current 
four-vectors, respectively. The dimensionless parameters T], A, and k charac- 
terize deviations from the standard model. (Since k gives a negligible con- 
tribution to double beta decay, we will not consider it from now on.) The 
coupling parameters i] and A, modified by the neutrino mixing, and denoted 
then usually as (77) and (A) are unknown (and presumably small). 

The lepton sector of the theory contains in general n generations of 
charged leptons as well as n left- and n right-handed neutrinos. The neu- 
trino mass matrix is the 2n x 2n matrix M 



6 Petr Vogel 



where Mo is the n x n lepton number conserving Dirac mass term, and 
the symmetric n x n matrices Ml and Mr are the lepton number violating 
Majorana mass terms. The matrix M has 2n real, but not necessarily positive, 
eigenvalues. Writing the eigenvalues as TUjCj, we can impose the physically 
reasonable condition that nij > 0. The sign of the eigenvalues of the mass 
matrix is contained in the phases ej = ±1 which are the intrinsic CP parities 
of the neutrinos j. 

Neutrino oscillation phenomena arise because the "mass eigenstates" of 
M or, more precisely their chiral projections and N^, are not necessarily 
the familiar weak interaction neutrinos that couple to the known intermedi- 
ate vector boson Wl and to the hypothetical right-handed boson Wr. The 
physical "weak eigenstate" or current neutrinos, the n left-handed neutrinos 
vl and the n right-handed ones v R (the prime has been added in order to 
stress that they are different particles), are related to the neutrinos of definite 
mass by the n x 2n mixing matrices U and V 



UN L , u' R = VN R . (9) 



The mixing matrices U and V obey the normalization and orthogonality 
conditions 

In 2n In 

YvtfJvi = 6 U , , Y,V$V Vi = S W , Y^U^Vvj = 0. (10) 
i=i i =1 i =1 

In neutrinoless (3j3 decay the rate depends on the effective parameters 
which are expressed in terms of the mixing matrices U and V: 

(A) = a ^ ' (; r,. ( r,. ( , (ii) 

3 
3 

(9v,x) ^X (9i,3 e i+ 9j,i e j) U e ,iUe,j ■ 

Here the prime indicates that the summation is over only relatively light 
neutrinos. Also, A and r\ are the dimensionless coupling constants for the 
right-handed current weak interaction, Eq.(^), and g^j are the coupling con- 
stants of interaction between the Majoron \ and the Majorana neutrinos Ni 
and Nj. For the heavy neutrino one obtains 

(m-^H = J2"^ m J 1U h . ( 12 ) 
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where the double prime indicates that the summation, involving the in- 
verse neutrino masses ttt," 1 , is over only the heavy neutrino mass eigenstates 
(rrij >1 GeV). 

It is now clear that, within the mechanism considered so far, there is no 
neutrinoless double beta decay if all neutrinos are massless. Not only (m„) 
vanishes in such a case but also (A) and (77) vanish due to the orthogonality 
condition Eq. ([h]). Moreover, (A) and (77) vanish for the same reason even if 
some or all neutrinos are massive but light and therefore the summation in 



Eq.(ll) contains all neutrino mass eigenstates. In that case, however, there is 
a smaller next order contribution from the mass dependence of the neutrino 
propagator, which for this purpose can be written as 



q 2 + m 2 




(13) 



The expression for e.g., (A) now contains J2j e jU e .jV e jm 2 which clearly 
shows that a nonvanishing neutrino mass is required. 

The presence of the phases tj in the expression for (m v ) means that 
cancellations are possible. In particular, for every Dirac neutrino there is an 
exact cancellation, since the Dirac neutrino is equivalent to a pair of Majorana 
neutrinos with the opposite sign of the phases €j and degenerate masses. 

In the general case the neutrinoless double beta decay rate is a quadratic 
polynomial in the unknown parameters 



m„r „ ,., (777.7. ) , „ , , (m„ 



[3%(0 + -> 0+)]- 1 = Ci±—£- + C 2 (A)-^-— ^ cos-01 + C 3 (77)-^—^ cos ^2 



(v —> u )\ — - r <^2W 

+ C 4 (X) 2 + C 5 ( V ) 2 + C 6 (X)(r,) cos(Vi - tA 2 ) . (14) 



Here ip\ and 7/^2 are the phase angles between the generally complex numbers 
m„, A and 77. (However, when CP invariance is assumed ^1,2 are either or 
7r.) The phase space integrals and the nuclear matrix elements are combined 



in the factors C 2 :. Assuming that we can calculate them, Eq. (14) represents 
an ellipsoid which restricts the allowed range of the unknown parameters 
(m„), (A) and (77) for a given value (or limit) of the OV double beta decay 
lifetime. 

In order to evaluate the nuclear matrix elements, we must consider the 
neutrino propagator. Assuming that (m u ) 2 is the only relevant quantity, one 
can perform the integration over the four-momentum of the exchanged par- 
ticle and obtain the "neutrino potential", which for m v < 10 MeV has the 
form 

where AE = (En) — l/2(Mj + Mf) is the average excitation energy of the 
intermediate odd-odd nucleus and the factor R (the nuclear radius) has been 
added to make the neutrino potential dimensionless. 
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When the Ov decay is mediated by the right-handed weak current in- 
teraction the evaluation of the decay rate becomes more complicated, since 
many more terms must be included (see |^,^,^|). If the four- momentum of 
the virtual neutrino is = to, q, the neutrino propagator contains 

w7o - q ■ 7 + rrij . 

The part of the propagator proportional to rrij is responsible for the neutrino 
potential Eq. (fl5l) . The part containing q leads to a new potential related to 
the derivative of H(r, AE), and the part with u> leads to yet another potential, 
which is a combination of H{r, AE) and its derivative. 

Similarly, there are now also more nuclear matrix elements, which contain 
in addition the nucleon momenta (i.e., the gradient operators), and depend 
on the nucleon spins and radii in a more complicated way (e.g., they contain 
tensor operators). The outgoing electrons are no longer just in the Sx/2 states, 
because for some of the operators one of the electrons will be in the pi/2 
state. The recoil matrix element, which originates from the recoil term in the 
nuclear vector current is numerically relatively large pp[ , resulting in more 
sensitivity to the parameter (77). The current best limits on (17) and (A) are 
listed in Q. 

4 Experimental techniques and results 

It is beyond the scope of this review to describe in detail the experimental 
techniques developed to meet the challenge of background suppression and 
signal recognition needed to determine the rate (or an interesting limit) of 
the (3(3 decay. Thus, only the briefest outline is given, and the most important 
experimental results are summarized in tables. 

Historically, the existence of f3f3 decay was first established using the geo- 
chemical method. Here one takes advantage of geologic integration times by 
searching for daughter products accumulated in ancient minerals that are 
rich in the parent isotope. (The related radiochemical method is applicable 
if the daughter isotope is radioactive.) Since the energy information is long 
lost, the mode of (3(3 decay responsible is not directly determined. Instead, 
the total decay rate is determined, and thus an upper limit of each mode as 
well. 

Only by measuring the energies of electrons released in the decay in the 
direct counting experiments can one distinguish directly the mode of decay. 
The 2v and Oi>x decay modes each result in a rather generic looking electron 
spectrum (see Fig.Q), and the observation of these decays requires either an 
extremely efficient background suppression or additional information, such 
as a tracking capability. 

The measured half-lives of the 2v mode are collected in Table Many of 
them have been measured by several groups; only the results with the smallest 
claimed errors are shown. (The case of 130 Te where the two competing results 
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have the same error but exclude each other is the only exception.) Also, 
the numerous half-life limits have been omitted. The 2i> mode is now well 
established; no doubt many more and more accurate results will become 
available soon. 

In fact, (3(3 decay is becoming a valuable tool of nuclear spectroscopy. 
The decay of 100 Mo into the excited 0+ state at 1130 keV in 100 Ru has been 
observed [ |30pl| ]. The technique used, the observation of the subsequent 7 
decay cascade, can be readily applied to other nuclei as well. This develop- 
ment not only expands the scope of the experimental study of the (3(3 decay, 
but allows more detailed comparison between theory and experiment (for an 
early attempt, see p^| ). 

Table 1. Recent (3(32 V results. 

(Only positive results are listed. The most accurate published values are given, 
except for 130 Te where two conflicting results with the same claimed errors are 
quoted.) 



Isotope T 2 / 2 (y) Reference 

[22] 
[23] 
[24] 
[25] 
[26] 
[27] 
1 27 1 
|2S| 
[25] 
|2S] 



geoch g eo chemical determination; total decay rate. 
radioch radiochemical determination; total decay rate 

The Ov mode can be approached quite differently from 2v and Ov, x modes 
because of the distinctive character of the fV electron sum spectrum — a 
monoenergetic line at the full Q- value (see Fig. |l|). Obviously, sharp energy 
resolution of 0V detectors is a big advantage which helps to isolate the line 
from background. As in the case of the 2v decay other features, such as 
tracking, naturally help as well. 

The best reported limits for the neutrinoless (3(3 decay modes are col- 
lected in Tables || and |[ Again, only the most restrictive limits for the given 
transition are shown. The longest half-life limit, reported for 76 Ge by the 



Ca (4.3±H ± 1.4) x 10 19 

Ge (1.77 ± 0.01 ±g i?) x 10 21 

Se (8.3 ± 1.0 ±0.7) x 10 19 

Zr (3.9 ± 0.9) x 10 19 geoch 

°Mo (6.82lo|3 ± °- 68 ) x 1018 

6 Cd (3.75 ±0.35 ±0.21) x 10 

8 Te (7.2 ± 0.4) x 10 24 3eoch 

°Te (2.7 ±0.1) x 10 21 3eoch 
(7.9 ± 1.0) x 10 20 3eoch 

°Nd (6.75±° f 2 ± 0.68) x 10 18 

8 U (2.0 ± 0.6) x 10 21 radioch 
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Heidelberg-Moscow collaboration [Q, is based on 24.16 kg yr of exposure 
and uses pulse shape discrimination to suppress the background (in the rele- 
vant energy region the background is a mere (0.06±0.02) events/(kg-yr-keV)). 
In that experiment 7 events were observed in the 3ct region around the 0v de- 
cay Q value, while from the background extrapolation one expects 13 events. 
Using this lack of background events an even more stringent limit (the entry 
in parenthesis in Table |^) is obtained. 

The limit based on the Te lifetime ratio in Table || is based on the different 
Q value dependence of the 0v and 2v modes. That this offers a valuable tool 
has been recognized already in the prophetic early paper by Pontecorvo |H . 
Even though the corresponding NME are not exactly equal, they are close 
enough to allow one to use the geochemical lifetime determination here and 
in Table |. 

Table 2. Best reported limits on Tj 1 ^ and (m„). 

The experimental result is listed first with its reference. This is followed by the limit 
on (m„) followed by the reference to the employed nuclear matrix element (NME). 
Whenever possible the choice of the authors of the experimental paper regarding 
the NME is respected. See the text for the discussion of uncertainties associated 
with the evaluation of NME's. 



Isotope 



T?/ 2 (10 22 y) (CL%) exp. ref. (m„) (eV) NME ref. 




is 



Ca 



7(i 



> 0.95 (76) 

> 1600(5700) (see 1} (90) 

> 2.7 (68) 

> 5.2 (68) 
Cd > 2.9 (90) 

T l/2(130) (see 2) 



; Ge 
2 Se 



Mo 



100 
116 



T 1/2 (128) 

136 Xe 

150 



°Nd 



(3.52 ±0.11) x 10 -4 

> 44 (90) 

> 0.12 (90) 



1 the first entry is based on the average background and the entry in parenthesis 
is based on the apparent lack of background counts in the corresponding energy 

interval. 

2 geochemical determination of the lifetime ratio 



5 Nuclear structure aspects 

The rate of the 2v (3(3 decay is simply 

l/T?Y 2 = G 2 „(E ,Z)\M 2v \ 2 , (16) 
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Table 3. The most restrictive Majoron limits. 



Isotope T "'* (y) and (CL %) (g UiX ) Reference 



48 Ca > 7.2 x 10 20 (90) < 5.3 x 10" 4 [41] 

76 Ge > 1.66 x 10 22 (90) < 1.8 x 10" 4 , [42] 

82 Se > 2.4 x 10 21 (68) < 2.3 x 10~ 4 [23| 

100 Mo > 5.4 x 10 21 (68) < 7.3 x 10" 5 [37| 

116 Cd > 1.2 x 10 21 (90) < 2.1 x 10 -4 [43] 

128 Te > 7.7 x 10 24 3eoch (90) < 3 x 10" 5 [27| 

136 Xe > 7.2 x 10 21 (90) < 1.6 x 10~ 4 [40] 

150 Nd > 2.8 x 10 20 (90) < 1 x 10~ 4 [ 25 1 



geoch 



geochemical determination; from total decay rate 



while for the neutrmoless decay (assuming that it is mediated by a light 
Majorana neutrino and that there are no right-handed weak interactions), 
and for the decay with Majoron emission, it is given by 

l/lf/a - GoAEo, Z)\M 0v \ 2 (m„) 2 , (17) 
1/T°;,f = G 0v , x (E ,Z)\M „, x \ 2 (g„, x ) 2 . 

Here the phase space functions G(Eq,Z) are accurately calculable, and the 
nuclear matrix elements M are the topic of this section. Obviously, the accu- 
racy with which the fundamental particle physics parameters (m„) and {g v ,x) 
can be determined is limited by our ability to evaluate these nuclear matrix 
elements. 

In that context there are three distinct set of problems: 

• 2v decay: the physics of the Gamow- Teller amplitudes 

• 0V decay with the exchange of light massive Majorana neutrinos: no selec- 
tion rules on multipoles, role of nucleon correlations, sensitivity to nuclear 
models. 

• 0V decay with the exchange of heavy neutrinos: physics of the nucleon- 
nucleon states at short distances. 



5.1 Two neutrino decay 

Since the energies involved are modest, the allowed approximation should 
be applicable, and the rate is governed by the double Gamow- Teller matrix 
element 

M iv = \p (/IK+IH x Hk r +IN) i / lg s 

GT ^ E m -(M i + M f )/2 [ > 
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where i, f are the ground states in the initial and final nuclei, and m are 
the intermediate 1 + (virtual) states in the odd-odd nucleus. The first factor 
in the numerator above represents the /3 + (or (n,p)) amplitude for the final 
nucleus, while the second one represents the f3~ (or (p, n)) amplitude for the 
initial nucleus. Thus, in order to correctly evaluate the 2v decay rate, we 
have to know, at least in principle, all GT amplitudes for both f3~ and /3 + 
processes, including their signs. The difficulty is that the 2v matrix element 
exhausts a very small fraction (1CP 5 — 10~ 7 ) of the double GT sum rule 0], 
and hence it is sensitive to details of nuclear structure. 

Various approaches used in the evaluation of the 2v decay rate have been 
reviewed recently in Ref. |6). The Quasiparticle Random Phase Approxima- 
tion (QRPA) has been the most popular theoretical tool in the recent past. 
Its main ingredients, the repulsive particle-hole spin-isospin interaction, and 
the attractive particle-particle interaction, clearly play a decisive role in the 
concentration of the (3~ strength in the giant GT resonance, and the rela- 
tive suppression of the (3 + strength and its concentration at low excitation 
energies. Together, these two ingredients are able to explain the suppression 
of the 2v matrix element when expressed in terms of the corresponding sum 
rule. 

Yet, the QRPA is often criticized. Two "undesirable" , and to some extent 
unrelated, features are usually quoted. One is the extreme sensitivity of the 
decay rate to the strength of the particle-particle force (often denoted as 
g pp ). This decreases the predictive power of the method. The other one is 
the fact that for a realistic value of g pp the QRPA solutions are close to 
their critical value (so called collapse). This indicates a phase transition, i.e., 
a rearrangement of the nuclear ground state. QRPA is meant to describe 
small deviations from the unperturbed ground state, and thus is not fully 
applicable near the point of collapse. Numerous approaches have been made 
to extend the range of validity of QRPA, see e.g. ||. Altogether, QRPA and 
its various extensions, with their ability to adjust at least one free parameter, 
are typically able to explain the observed 2v decay rates. 

At the same time, detailed calculations show that the sum over the excited 
states in Eq.(|l8|) converges quite rapidly In fact, a few low lying states 
usually exhaust the whole matrix element. Thus, it is not really necessary to 
describe all GT amplitudes; it is enough to describe correctly the f3 + and (3~ 
amplitudes of the low- lying states, and include everything else in the overall 
renormalization (quenching) of the GT strength. 

Nuclear shell model methods are presently capable of handling much 
larger configuration spaces than even a few years ago. Thus, for many nuclei 
the evaluation of the 2v rates within the Qfiui shell model space is feasi- 
ble. (Heavy nuclei with permanent deformation, like 150 Nd and 238 U remain, 
however, beyond reach of the shell model techniques.) Using the shell model 
avoids, naturally, the above difficulties of QRPA. At the same time, the shell 
model can describe, using the same method and the same residual interaction, 



Double beta decay 13 



a wealth of spectroscopic data, allowing much better tests of its predictive 
power. 

5.2 Neutrinoless decay: light Majorana neutrino 

If one assumes that the Ov decay is caused by the exchange of a virtual light 
Majorana neutrino between the two nucleons, then several new features arise: 
a) the exchanged neutrino has a momentum q ~ l/r nn ~ 50 — 100 MeV (r nn 
is the distance between the decaying nucleons). Hence, the dependence on 
the energy in the intermediate state is weak and the closure approximation is 
applicable and one does not have to sum explicitly over the nuclear interme- 
diate states. Also, b) since qR > 1 (R is the nuclear radius), the expansion 
in multipolcs is not convergent, unlike in the 2i> decay. In fact, all possi- 
ble multipoles contribute by a comparable amount. Finally, c) the neutrino 
propagator results in a neutrino potential of a relatively long range (see Eq. 
(15)). 

Thus, in order to evaluate the rate of the 0v decay, we need to evalu- 
ate only the matrix element connecting the ground states + of the initial 
and final nuclei. Again, we can use the QRPA or the shell model. Both cal- 
culations show that the features enumerated above are indeed present. In 
addition, the QRPA typically shows less extreme dependence on the particle- 
particle coupling constant g pp than for the 2v decay, since the contribution 
of the 1 + multipole is relatively small. The calculations also suggest that for 
quantitatively correct results one has to treat the short range nucleon-nucleon 
repulsion carefully, despite the long range of the neutrino potential. 

Does that mean that the calculated matrix elements are insensitive to nu- 
clear structure? An answer to that question has obviously great importance, 
since unlike the 2v decay, we cannot directly test whether the calculation is 
correct or not. 

For simplicity, let us assume that the Qv (3/3 decay is mediated only by the 
exchange of a light Majorana neutrino. The relevant nuclear matrix element 
is then the combination Mq T — M F V , where the GT and F operators change 
two neutrons into two protons, and contain the corresponding operator plus 
the neutrino potential. One can express these matrix elements either in terms 
of the proton particle - neutron hole multipoles (i.e., the usual beta decay 
operators) or in terms of the multipole coupling of the exchanged pair, nn 
and pp. 

When using the decomposition in the proton particle - neutron hole mul- 
tipoles, one finds that all possible multipoles (given the one-nucleon states 
near the Fermi level) contribute, and the contributions have typically equal 
signs. Hence, there does not seem to be much cancellation. 

However, perhaps more physical is the decomposition into the exchanged 
pair multipoles. There one finds, first of all, that only natural parity multi- 
poles (it = (— l) 7 ) contribute noticeably. And there is a rather severe cancel- 
lation. The biggest contribution comes from the + multipole, i.e., the pairing 
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Fig. 2. The cumulative contribution, i.e., the summed contribution of all natural 
parity multipoles up to I of the exchanged nn and pp pair, to the Ov nuclear matrix 
element combination M^ T — Mp v . The full line is for 76 Ge and the dashed line for 



part. All other multipoles, related to higher seniority states, contribute with 
an opposite sign. The final matrix element is then a difference of the pairing 
and higher multipole (or broken pair = higher seniority) parts, and is consid- 
erably smaller than either of them. This is illustrated in Fig. ^| where the cu- 
mulative effect is shown, i.e., the quantity M(I) = J2j [ m gt( j ) ~ m f"( j )] 
is displayed for 76 Ge (from j47|) and 48 Ca (from Q). Thus, the final result 
depends sensitively on both the correct description of the pairing and on the 
admixtures of higher seniority configurations in the corresponding initial and 
final nuclei, ft appears, moreover, that the final result might depend on the 
size of the single particle space included. That important question requires 
further study. 

Since there is no objective way to judge which calculation is correct, one 
often uses the spread between the calculated values as a measure of the 
theoretical uncertainty. This is illustrated in Fig. ||[ There, I have chosen two 
representative QRPA sets of results, the highly truncated "classical" shell 
model result of Haxton and Stephenson, and the result of more recent shell 
model calculation which is convergent for the set of single particle states 
chosen (essentially Ofrui space). 

For the most important case of 76 Ge the calculated rates differ by a factor 
of 6-7. Since the effective neutrino mass (m„) is inversely proportional to the 
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Fig. 3. Half-lives (in years) calculated for (m„) = 1 eV by various representative 
methods and different authors for the most popular double-beta decay candidate 
nuclei. Solid lines are QRPA from |$5|, dashed lines are QRPA from j$9| (recalcu- 
lated for qa = 1-25 and a' = —390 MeV fm 3 ), dotted lines are shell model and 
dot-and-dashed lines are shell model 1491. 



square root of the lifetime, the experimental limit of 1.6 x 10 25 y translates 
into limits of about 1 eV using the NME of |ff and about °- 4 eV with 
the NME of [p|j35|. On the other hand, if one would accept the more stringent 
limit of 5.7 x 10 25 [^3), even the more pessimistic matrix elements restrict 
(m„) < 0.5 eV. Needless to say, a more objective measure of the theoretical 
uncertainty would be highly desirable. 

In Tables | and | we list the deduced limits on the fundamental param- 
eters, the effective neutrino Majorana mass (m v ), and the Majoron coupling 
constant {g v ,x)- The references to the source of the corresponding nuclear 
matrix elements, used to translate the experimental half-life limit into the 
listed limits on (m„) and {g v ,x) are a l so given. When using the tables one 
has to keep in mind the uncertainties illustrated in Fig. El. 

5.3 Neutrinoless decay: very heavy Majorana neutrino 

The neutrinoless (3(3 decay can be also mediated by the exchange of a heavy 
neutrino. The decay rate is then inversely proportional to the square of the 
effective neutrino mass [^o| . In this context it is particularly interesting to 
consider the left-right symmetric model proposed by Mohapatra [^l|. In it, 
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Fig. 4. The Feynman graph description of the parity-violating nncleon-nucleon 
force (left graph) and of the (3/3 decay with the exchange of a heavy neutrino 
mediated by the pion exchange. The short range lepton number violating amplitude 
is symbolically described by the filled blob in the right graph. 



one can find a relation between the mass of the heavy neutrino Mn and the 
mass of the right-handed vector boson Wr. Thus, the limit on the (3(5 rate 
provides, within that specific model, a stringent lower limit on the mass of 
W R . 

The process then involves the emission of the heavy by the first neu- 
tron and its virtual decay into an electron and the heavy Majorana neutrino, 
— > e~ + vn- This is followed by the transition vn — > e~ + W£ and the 
absorption of the on the second neutron, changing it into the second pro- 
ton. Since all exchanged particles between the two neutrons are very heavy, 
the corresponding "neutrino potential" is of essentially zero range. Hence, 
when calculating the nuclear matrix element, one has to take into account 
carefully the short range nucleon-nucleon repulsion. 

As long as we treat the nucleus as an ensemble of nucleons only, the only 
way to have nonvanishing nuclear matrix elements for the above process is 
to treat the nucleons as finite size particles. In fact, that is the standard way 
to approach the problem f30f ; the nucleon size is described by a dipole form 
factor with the cut-off parameter A ~ 0.85 GeV. Using such a treatment of 
the nucleon size, and the half- life limit for the 76 Ge 0v decay listed in Table 
pj, o ne obtains a very interesting limit on the mass of the vector boson Wr 

mw R > 1.6 TeV . (19) 

However, another way of treating the problem is possible, and already 
mentioned in |5(| . Let us recall how the analogous situation is treated in the 
description of the parity- violating nucleon-nucleon force pa. There, instead 
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of the weak (i.e., very short range) interaction of two nucleons, one assumes 
that a meson (jr,ui,p) is emitted by one nucleon and absorbed by another 
one. One of the vertices is the parity-violating one, and the other one is the 
usual parity-conserving strong one. The corresponding range is then just the 
meson exchange range, easily treated. The situation is schematically depicted 
in the left-hand panel of Fig. |5.3| . The analogy for [3(3 decay is shown in 
the right-hand graph. It involves two pions, and the "elementary" lepton 
number violating 0j3 decay then involves a transformation of two pions into 
two electrons. Again, the range is just the pion exchange range. It would be 
interesting to see if a detailed treatment of this graph would lead to more 
or less stringent limit on the mass of the Wr than the treatment with form 
factors. The relation to the claim in (Hi) that an analogous graph contributing 
to the lepton number violating muon capture identically vanishes should be 
further investigated; in fact that claim is probably not valid. 

6 Perspectives 

As stated earlier, the present best limits on the rate of the Ov (3f3 decay, 
or equivalently on the neutrino effective Majorana mass (m u ), were obtained 
with an exposure of about 20 kg-yr. Several experiments (Heidelberg-Moscow 
76 Ge ||, IGEX 76Ge II' Caltech-Neuchatel TPC 136 Xe @) are presently 
at or near that level. The other limits in Tables || and || were obtained with 
smaller exposures of ~ 1 kg-yr. The detector NEMO-3, with planned source 
mass of 10 kg, is being built and should be operational soon In few years 
of operation it should reach a half-life limit of ~ 10 25 yr for the Qv decay 
of 100 Mo, and perhaps other nuclei as well. However, further improvements 
with the existing detectors becomes increasingly difficult, since the sensitivity 
to the (m u ) is proportional to only the 1/4 power of the source mass and 
exposure time. Thus, much larger sources are clearly needed. 

What are the perspectives of a radical improvement in the search for the 
0v /?/? decay? To achieve that, one would have to build a detector capable of 
using hundreds of kg or even several tons of the source material. At the same 
time, the background per unit mass has to be correspondingly improved so 
that one can benefit from the larger mass. Obviously, such program is very 
challenging. 

The difficulty begins with the problem of acquiring such a large mass of the 
isotopically separated and radioactively clean material. Here, the principal 
obstacle is the cost of the isotope separation. (This can be avoided only if 
the source isotope has large abundance; in practice that is true only for 130 Te 
with 34% abundance.) 

The second unavoidable difficulty is the background caused by the 2v 
decay. One can observe the 0v decay only if its rate exceeds the fluctuations 
of the 2v events at the same energy, i.e., near the decay Q value. The number 
of 2v decays in an energy interval AE near Q depends on these quantities 
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like ~ (AE/Q) 6 provided AE <C Q. (If the energy resolution is folded in, 
this dependence is somewhat modified.) Thus, good energy resolution, which 
determines how wide interval AE one must consider, is again crucial in order 
to reduce the effect of this "ultimate" background. 

One of the proposal for such a large (3(3 experiments has been extensively 
discussed in the literature (see, e.g., 0). The project, with the acronym GE- 
NIUS, would use a large amount of 'naked' enriched 76 Ge, in the form of 
an array of about 300 detectors, suspended in liquid nitrogen, which pro- 
vides simultaneously cooling and shielding. It is envisioned that the detector 
would consist of one ton of enriched 76 Ge. The anticipated background is 
0.04 counts/ (keV-yr-t), i.e., about 1000 times lower than the best existing 
backgrounds. Such a detector could reach the half-life limit of about 6 x 10 27 
yr within one year of operation, thus improving the neutrino mass limit by 
an order of magnitude. 

Another large project, CUORE, |57| is a cryogenic set-up consisting of 17 
towers, each containing 60 cubic crystals of Tc02- It would be housed in a 
single specially constructed dilution refrigerator and would contain about 800 
kg of the sensitive material. A prototype system, CUORICINO, consisting of 
one of the towers, is being developed now. 

An experiment with a large amount (100 tons) of natural molybdenum, 
(abundance of the (3(3 candidate 100 Mo = 9.6%) with good energy and posi- 
tion resolution, is proposed in Ref. Q. 

In order to radically suppress the background, the Ba ions, the final prod- 
ucts of the 136 Xe double beta decay, could be identified by laser tagging. That 
approach, described in Ref. | p9[ , would allow to use a large Time Projection 
Chamber with perhaps ton quantities of 136 Xe, reaching sensitivities to half- 
lives ~ 10 28 years. 

This, still incomplete list of proposed very large (3(3 decay experiments 
shows that the field is entering a critical phase. If the new techniques, men- 
tioned above, can be developed in conjuntion with the large source mass, the 
background caused by radioactivities can be essentially eliminated. However, 
as stated above, the ultimate background due to the tail of the 2v decay can 
be compensated only by a superior energy resolution. 

Given the importance of the neutrinoless decay, it is likely that several 
of these large and costly projects, involving ton-years of exposure and a cor- 
respondingly reduced background, will be realized in the foreseeable future. 
Thus, sensitivity to the neutrino Majorana mass (m u ) approaching 0.01 eV 
may be in sight. Whether the neutrinoless decay will be discovered is un- 
known, but the reasons to look for it are so compelling, that the search will 
undoubtedly continue. 
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7 Implications 

The study of Ov (3(3 decay provides at present an upper limit well below 1 eV 
for the effective electron neutrino Majorana mass (m„) even if the most pes- 
simistic nuclear matrix elements are used. What are the consequences of that 
limit when combined with the manifestations of the neutrino oscillations? 

Recall that the atmospheric neutrino anomaly (with its zenith angle de- 
pendence) implies the nearly maximum mixing of (i and r neutrinos (or \x 
and sterile neutrinos) with Am 2 ~ 1CP 3 eV 2 (see chapter 5 of this book). 
There is, so far, no unique neutrino oscillation solution to the solar neutrino 
dehcit(see chapter 4 of this book). However, all of the acceptable solutions 
have Am 2 < 10 -4 eV 2 and involve electron neutrinos. Both large and small 
mixing angle solutions are currently compatible with the data. Finally, the 
third "positive" evidence comes from the LSND experiment (see chapter 7 
of this book), and implies relatively small mixing between the electron and 
muon neutrinos and Am 2 > 0.1 eV 2 . A full analysis must contain, in addition, 
all experimental results which exclude various parts of the possible regions of 
the quantities Am 2 and the mixing angles. Taking all these findings together 
would necessarily imply the existence of a fourth neutrino, which must be 
"sterile" given the constraint on the invisible width of the Z. At the same 
time, it is well known that oscillation experiments are not able to furnish the 
overall scale of the neutrino masses. 

This absolute neutrino mass scale is essential not only as a matter of 
principle, but in particular if one wants to ascribe part of the dark matter, 
namely its "hot" component, to massive neutrinos. Doing that would mean 
that the sum of the neutrino masses ^2 m v is one or several eV. Tritium beta 
decay gives an upper limit of a similar magnitude for any mass eigenstate with 
a large electron flavor component. Clearly, if light neutrinos are responsible 
for a nonnegligible part of the dark matter, the oscillation data mean that 
at least two, and possibly all neutrino masses are nearly degenerate. (Such 
scenario was discussed for the first time in Ref. pQ).) The relation of the Op 
(3(3 decay and the oscillation scenarios, in particular the scenarios involving 
degenerate neutrinos, has been a topic of several recent papers luTjo^Jr^J&if . 

The consequences are particularly dramatic if one assumes that only three 
massive Majorana neutrinos exist with nearly degenerate masses m, ~ m~ 
O(eV). (Hence discarding for this purpose the LSND experimental result, 
even though there is no evidence against it.) The (3(3 decay constraint can be 
expressed as Q 

(m„) = |m lC 2 c 2 e^ + m 2 c 2 s 2 e^' + m 3 s 2 e l25 | , (20) 

where c% , Sj denote cos 9i , sin 9i in the 3x3 mixing matrix, 5 is the CP violating 
phase in that matrix, and <f>' are the CP violating phases in the diagonal 
mass matrix. Clearly, the differences in m, can be neglected in this case. 
Moreover, the reactor long baseline experiments have established that v e do 
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not mix very much with anything else near Am 2 ~ 10~ 3 eV 2 , which means 
that the angle 62 is small. At the same time, the angle 9\ which controls 
the atmospheric neutrino oscillations is near its maximum value sin 2 29 \ ~ 1. 
Thus 

|cos 2 6»3 + sin 2 <9 3 e' iW '- 0) | < (m 1/ )/m<l . (21) 

Hence also #3 must be near the maximum mixing, sin 26% ~ 1, and the CP 
phases in the above equation (|2l| ) are such that the two terms cancel each 
other. 

That would be a very unexpected result. We would have three massive 
highly degenerate neutrinos with bimaximal mixing. Moreover, the electron 
neutrino would be 'quasi-Dirac' with its two components essentially canceling 
each other in their contribution to the Ov (3(3 decay. While such a scenario is 
rather problematic (see |52|), and it does not accommodate the LSND result 
at all, it illustrates the power of the neutrinoless (3(3 decay in constraining 
the choice of the neutrino oscillation scenarios. 

8 Conclusions 

The quest for neutrino mass is at a critical stage at present. The evidence for 
neutrino mixing is getting stronger and stronger, and the basic parameters 
describing the neutrino oscillation phenomena are being constrained more 
and more. At the same time, the oscillation searches cannot give us the scale 
of the neutrino masses, but only their differences. Among the experiments 
that are sensitive to the masses themselves, albeit to their different aspects 
(end point of the ordinary beta decay, observation of the supernova neutrinos, 
and the neutrinoless double beta decay), only the Ov decay is able to reach 
the sub eV region, and in a foreseeable future extend it by a substantial 
margin. 

In this review the present status of the (3(3 decay is described. The un- 
pleasant uncertainty related to the nuclear structure aspect of the problem is 
estimated to be at the level of a factor of 2-3 for the effective neutrino mass. 
However, the experimental progress is such that even using the most con- 
servative nuclear matrix elements allows us to push the limit well below the 
competing techniques. The nuclear structure uncertainty can be reduced by 
further development of the corresponding nuclear models. At the same time, 
by reaching comparable experimental limits in several nuclei, the chances of 
a severe error in the NME will be substantially reduced. 

Several projects are under way that will improve the life-time limit sub- 
stantially, or find the Ov decay. Already now the search for (3(3 decay gives 
important constraints on the fundamental properties of neutrinos and their 
interactions. The role of the (3(3 decay in the whole enterprise descibed in 
various chapters of this book will be substantially strengthened once these 
ambitious projects are underway. 
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